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Abstract. We present a set of detailed, self-consistent, isotropic dy- 
namical models for disc galaxies. We start from the hypothesis that 
each galaxy can be decomposed in a bulge, following the r 1 / 4 law, and a 
disc with an exponential projected density profile; and that the isoden- 
sity surfaces of each component can be represented by similar concentric 
spheroids. Under these conditions we produce the rotational velocity and 
velocity dispersion profiles, after taking into account both the asymmetric 
drift effects and the integration along the line of sight. 

The model is successfully applied to reproduce the stellar kinematic 
and photometry of the bulge of the 4 lenticular and early-type spiral 
galaxies NGC 4565, NGC 7814, NGC 5866 and NGC 4594. For these 
galaxies detailed stellar kinematical data are available at different posi- 
tions across the galaxy. 

The application of our models shows that: 1) For all the galaxies 
considered in this work, an isotropic model is able to reproduce the whole 
dynamical data. This is surprising in the special case of NGC 4565, for 
which previous models were unable to reproduce the velocity distribution 
without introducing strong anisotropies. 2) We do not need a dark mass 
halo to reproduce these data. This do not mean that a dark halo is not 
present, but just that its dynamical effects in these inner regions of the 
galaxies are negligible. 

In our opinion, these results strongly support our guess that a com- 
plete dynamical model, including both a disc and a bulge component, is 
needed to understand the structure of the lenticular and Sa galaxies. 



1. Introduction 

The study of the distribution function of stars in galaxies may give strong con- 
straints on the scenarios of formation of the galaxies. In the past, a particular 
attention has been paid to distinguish between isotropic cases, where the distri- 
bution function is characterized by only two classical integrals of motion (energy 
and angular momentum) and the anisotropic case, where a third (unknown) iso- 
lating integral is needed to explain the dynamics of the system. 

A widely used tool to investigate the distribution function in elliptical galax- 
ies is the classical V/a parameter, first introduced by Binney (1978). By ap- 
plying this method to a large number of elliptical galaxies, it seems that two, 
different kind of ellipticals may exist: elliptical galaxies fainter than M#=19.5, 
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Figure 1. Photometry of a model with two components having 
different axial ratios. We assume a spherical bulge and an exponential 
disk with b/a = 0.1. On the left figure we show the isophotes of such 
a model, while on the right is shown the behavior of the global axial 
ratio as a function of the radius. Notice how the isophotes becomes 
more and more flatter while moving out from the center, and then turn 
back to be round when the bulge component become dominant again. 
All radii are in units of R P . 



with outer isophotes quite flat (disky), which are represented well by isotropic 
models, and bigger ellipticals (Mb < 19.5), with isophotes more "boxy" and 
that cannot be explained by a model with isotropic velocity ellipsoid. 

The structure and isotropy of the bulges of disk galaxies is far less clear: 
working on the same set of data, concerning fours SO and Sa galaxies, Kormendy 
& Illingworth (1982) claims the isotropy of their bulges, while Whitmore, Rubin 
and Ford (1985) deduce the opposite conclusion. The fact that the same data 
may be interpreted in such different ways clearly demonstrates the ambiguity of 
the V/a test when applied to the bulges. 

On the other hand, the informations about the dynamical structures of 
the bulges could be crucial to determine if bulges can be considered as "dwarf 
ellipticals embedded in a gas disk" , as stellar populations, colors and luminosity 
profiles seems to suggest, or if they are a totally different class of objects. 



2. The model 

The main problems encountered dealing with models of bulges is the mixing 
of stars between disk and bulge and the mutual interaction of their potentials. 
These effects make harder the problem of to study the "pure" bulge component, 
that require some preliminary assumptions. 
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Figure 2. Rotational curves and velocity dispersion profiles for mod- 
els with different values of the parameters Rb/Rd, Mb/Mjj. Here we 
adopted the values Mb/Md = 1.36 which, according to Simien and 
de Vaucouleurs (1986), corresponds to a typical Sa galaxy, while the 
Rb/Rd parameter is ranging from 0.1 to 0.9 (models Sal through Sa9). 



In modeling the elliptical galaxies it is common to assume that the isodensity 
surfaces are similar, concentric spheroids. This is in reasonable agreement with 
the shape of the observed isophotes, and has the great advantage to reduce the 
complexity of the dynamical equations. As said before, in such a case (ellipsoidal 
structures) the V/a test can be applied. 

On the other hand, the isophotes of bulges in disk galaxies shows many 
different behaviors, such as wide change of flattening from the center of the 
galaxy to the outer regions. In this case, it is not possible to apply a simple 
ellipsoidal model. This problem is particularly presenting SO and Sa galaxies, 
where the contributions of the bulge- and disk- components are comparable in 
strength. 
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Starting from a previous work of Galletta et al. (1990), we extended the 
ellipsoidal hypothesis to early-type disk galaxies, assuming that each galaxy 
may be "splitted" in two different components, having each one the following 
properties: 

• The isodensity surfaces are similar concentric spheroids; 

• The anisotropy parameter j3 is constant throughout the galaxy (Here, we 
present only the case j3 = of the isotropic model) . 

• The galaxy is rotating about the z-axis with velocity V(R, Z), consistent 
with the self-gravity hypothesis. 

• Each component has a constant M/L ratio. 

Under these hypothesis, the circular velocity due to the potential of each 
component is given by the well-known form (see Binney & Tremaine, 1986): 

v 2 (R, z) = -2<kGVT^R P(m2)dT (1) 

cV ' Jo (r + l)Vr + 1 - e 2 V ' 



with e the ellipticity of the component, and m defined as 

2 R 2 Z 2 

m 2 = + ? (2) 



The density profile of the bulge is assumed to follow the r 1 / 4 law, while the 
disc (also "ellipsoidal", but much flatter than the bulge) is described with an 
exponential density profile. In Fig. |l] we present the photometric behavior of a 
typical galaxy constructed with these assumptions. 

Since we are dealing with the inner kinematical properties of the galaxies, 
we deliberately discard any dark matter contribution to the overall potential, to 
limit the number of free parameters. 



3. A library of models 

A first goal of this work was to find out how kinematical (rotation and velocity 
dispersion) curves are changing when different bulge/disk ratios are assumed. 
Under the hypotheses described in the previous section, we produced several 
different models for different values of the parameters. In Fig. Q and in Fig. ||] 
some of the obtained results are shown. 
We want to stress out some points: 

• The rotation curves of the galaxies on a plane perpendicular to the main 
galaxy plane become more and more linear at higher distance from this 
latter plane. This is in agreement with the available data for most of the 
observed spiral bulges; 

• On the other hand, the asymmetric drift effect is very strong in the nucleus; 
this is due to the steep increase of density for the inner regions of the bulge 
(described as a pure r 1 / 4 model, without any core radius). 
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• The velocity dispersion profile is always not gaussian; more likely, it can 
be approximated by the sum of two gaussian curves. 
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Figure 3. Behavior of the rotational velocity and velocity dispersion 
profile in a pure bulge model with a given axial ratio b/a = 0.6. The 
different profiles are extracted at different vertical shifts z from the 
major axis, as shown in the legend. Both the asymmetric drift effect 
and the line-of-sight projection are taken into account. All the models 
are seen edge-on and along the major axis. 



4. Comparison to the observations 

As a second goal, we applied our model to the galaxies studied by Kormendy 
and Illingworth (1982) and by Whitmore et al. (1985). In Tab. 1 the best-fit 
parameters found for all the four galaxies are shown. Our isotropic model fit 
fairly well the observations; this is particularly amazing in the case of NGC 4565 
(Fig. |j), for which both the papers conclude that strong anisotropies in the 
velocity distribution are present. 
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Figure 4. The rotational curves and velocity dispersion profiles 
predicted by the model compared with the observations of Kormendy 
& Illingworth (1982). The label || means that the data are extracted 
along an axis which is parallel to the major axis, while the label _L 
label means parallel to the minor axis. All the velocities and velocity 
dispersions are in Km/sec. As described into the text, the model is 
described by one free parameter only. 

The advantage of our model is that the only free parameter left is the M/L 
ratio of the stellar component, if a detailed photometry is available. In such 
a case, the decomposition of the galaxy into a bulge and a disk component 
can be properly made, and we are able to predict the observed kinematics and 
photometry at any point of the galaxy on the sky simply by changing the value 
of this parameter. The diagnostic power of the model is, under these conditions, 
very high. 
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Table 1. Best-fit parameters of the model for the four early- type disk 
galaxies studied by Kormendy & Illingworth. 





NGC 4565 


NGC 4594 


NGC 5866 


NGC 7814 


Hubble Type 


Sb 


Sa 


SO 


Sab 


Lb/Ld 


0.538 


11.5 


9 


1.63 


Rb 


50.59" 


18.9" 


32.2" 


31.3" 


Rd 


122.5" 


90" 


41" 


117" 


(b/a) B i 


0.59 


0.55 


0.47 


0.53 


D(Mpc) 


18.9 


19.3 


18.4 


25.0 


M B 


-21.11 


-22.26 


-20.47 


-20.73 


(M/L) b 


9.5 


7.5 


7. 


12. 


M tot 


4.17 


8.6 


1.72 


4.04 


(b/a) B 


0.59 


0.55 


0.8 


0.55 


Rb/Rd 


0.7 


1.0 


1.0 


0.487 


M B /M D 


1.06 


1.0 


9.0 


1.63 



5. Conclusions 

We developed a dynamical model to study the presence of isotropy in the bulges 
of early-type disk galaxies. The main advantage on previous works is that both 
the bulge and disk contributions may be properly taken into account, as well as 
the presence of the asymmetric drift and the integration along the line-of-sight. 
The application of our model shows that: 

• For all the galaxies considered in this work, an isotropic model is able to 
reproduce the whole dynamical data. This is surprising in the special case 
of NGC 4565, for which previous works suggested strong anisotropics in 
the velocity distribution. 

• The model don't needs a dark halo to reproduce these data. Its presence 
is, of course, possible, but its dynamical effects are negligible in the inner 
regions ( R < 1.2R e ) of the considered galaxies. This result is in agreement 
with the conclusions of other authors (e.g. Persic & Salucci, 1995). 

• The method seems to have a great diagnostic power if good kinematical 
data at different position angles are available and/or if many offset spectra 
have been taken. 
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